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three-dimensional Euler equations 1 
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Abstract 

More than 160 years after their invention by Hamilton, quaternions are now widely used 
in the aerospace and computer animation industries to track the orientation and paths of 
moving objects undergoing three-axis rotations. It is shown here that they provide a natural 
£S| \ way of selecting an appropriate ortho-normal frame - designated the quaternion-frame - for 

5^ ' a particle in a Lagrangian flow, and of obtaining the equations for its dynamics. How these 

ideas can be applied to the three-dimensional Euler fluid equations is then considered. This 
work has some bearing on the issue of whether the Euler equations develop a singularity 
■ in a finite time. Some of the literature on this topic is reviewed, which includes both the 

Beale-Kato-Majda theorem and associated work on the direction of vorticity by Constantin, 
Fefferman & Majda and Deng, Hou and Yu. It is then shown how the quaternion formulation 
■^j- ■ provides an alternative formulation in terms of the Hessian of the pressure. 

O ' 

o 
o 

\q ' This paper is dedicated to the memory of Victor Yudovich (1934-2006) with whom the author 
O ■ discussed some of these ideas in their early stages. 
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1 General introduction 

1.1 Historical remarks 

Everyone loves a good story : William Rowan Hamilton's feverish excitement at the discovery 
of his famous formula for quaternions on 16th October 1843 as a composition rule for orienting 
his telescope; his inscription of this formula on Broome (Brougham) Bridge in Dublin; and then 
his long and eventually unfruitful championing of the role of quaternions in mechanics, are all 
elements of a story that has lost none of its appeal [1, 2]. Hamilton's name is still revered today 
for the audacity and depth of his ideas in modern mechanics and what we now call symplectic 
geometry [3, 4, 5]. Indeed, evidence of his thinking is everywhere in both classical and quantum 
mathematical physics and applied mathematics, yet in his own century his work on quater- 
nions evoked criticism and even derision from many influential fellow scientists 3 . Ultimately 
quaternions lost out to the tensor notation of Gibbs, which is the basis of the 3- vector notation 
universally used today. 

In essence, Hamilton's multiplication rule for quaternions represents compositions of rota- 
tions [6, 7, 8, 9, 10, 11]. This property has been ably exploited in modern inertial guidance 
systems in the aerospace industry where computing the orientation and the paths of rapidly 
moving rotating satellites and aircraft is essential. Kuipers' book explains the details of how 
calculations with quaternions in this field are performed in practice [12]. Just as importantly, 
the computer graphics community also uses them to determine the orientation of tumbling ob- 
jects in animations. In his valuable and eminently readable book, Andrew Hanson [2] says in 
his introduction : 

Although the advantages of the quaternion forms for the basic equations of attitude control - 
clearly presented in Caylcy [6] , Hamilton [7, 8] and especially Tait [9] - had been noticed by 
the aeronautics and astronautics community, the technology did not penetrate the computer 
animation community until the land-mark Siggraph 1985 paper of Shoemake [13]. The 
importance of Shoemake's paper is that it took the concept of the orientation frame for 
moving 3D objects and cameras, which require precise orientation specification, exposed 
the deficiencies of the then-standard Euler-angle methods , and introduced quaternions to 
animators as a solution. 

Hamilton's 19th century critics were, of course, correct in their assertion that quaternions need 3- 
vector algebra to manipulate them, yet the use the aero/astronautics and animation communities 
have made of them are one more illustration of the universally acknowledged truth that while 
new mathematical tools may not be of immediate use, and may appear to be too abstract or 
overly elaborate, they may nevertheless turn out to have powerful applications undreamed of at 
the time of their invention. 

1.2 Application to fluid dynamics 

The close association of quaternions with rigid body rotations [9, 10, 11] points to their use 
in the incompressible Euler equations for an inviscid fluid as a natural language for describing 
the alignment of vorticity with the eigenvectors of the strain rate that are responsible for its 

3 Kelvin was one such example: see [1]. 

4 A well-known deficiency of Euler-angle methods lies in the problems they suffer at the poles of the sphere 
where the azimuthal angle is not defined. 
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nonlinear evolution. For a three-dimensional fluid velocity field u(x, t) with pressure p(x, t), 
the incompressible Euler equations are [14, 15, 16, 17, 18] 

Du 

jh—Vp, (1.1) 
where the material derivative is defined by 

The motion is constrained by the incompressibility condition div u = 0. The crucial dynamics lies 
in the evolution of the velocity gradient matrix Vit = {uij} which comes from the differentiation 
of (1-1) ^ 

—jjf = ~ u i,kUk,j ~ Pij , (1-3) 

where Py is the Hessian matrix of the pressure 

P* = ■ (1-4) 

OXiOXj 

The incompressibility condition divw = insists that Truij = which, when applied to (1.3), 
gives 

TrP = Ap = -u ijk u kji = 5 
In (1.5) above, S is the strain matrix whose elements are defined by 



Tr P = Ap = -u iik u k!i = \u? - Tr (S 2 ) . (1.5) 



S ij = I ( u i,J + u J,i) ■ ( L6 ) 

This is a symmetric matrix the alignment of whose eigenvectors e, is fundamental to the dynam- 
ics of the Euler equations. For instance, vortex tubes and sheets (Burgers' vortices and shear 
layers) always have one eigenvector aligned with the vorticity vector u [18]. 

This review cannot hope to deal with every aspect of the three-dimensional Euler equations, 
particularly the vast literature on weak and distributional solutions; the reader is urged to read 
the book by Majda & Bertozzi [15] to study these aspects of the problem. Here we concentrate 
on one particular aspect, which is the role played by quaternions in providing a natural language 
for extracting geometric information from the evolution of Uij. Because they are particularly 
effective in computing the orientation of rotating objects moving in three-dimensional paths they 
might be useful in understanding how general Lagrangian flows behave, particularly in finding 
the evolution of the ortho-normal frame of particles moving in such a flow. These particles could 
be of the passive tracer type transported by a background flow or they could be Lagrangian fluid 
parcels. Recent experiments in turbulent flows can now detect the trajectories of tracer particles 
at high Reynolds numbers [19, 20, 21, 22, 23, 24, 25, 26, 27, 28] : see Figure 1 in [19]. For any 
system involving a path represented as a three-dimensional space-curve, the usual practice is 
to consider the Frenet-frame of a trajectory constituted by the unit tangent vector, the normal 
and the bi-normal [2, 28]. In navigational language, this represents the corkscrew- like pitch, 
yaw and roll of the motion. While the Frenet-frame describes the path, it ignores the dynamics 
that generates the motion. Attempts have been made in this direction using the eigenvectors 
ej of S but ran into difficulties because the equations of motion for are unknown [29]. In §2 
another ortho-normal frame is introduced that is associated with the motion of each Lagrangian 
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particle. It is designated the quaternion-frame : this frame may be envisioned as moving with 
the Lagrangian particles, but its evolution derives from the Eulerian equations of motion. The 
advantage of this approach lies in the fact that the Lagrangian dynamics of the quaternion-frame 
can be connected to the fluid motion through the pressure Hessian P defined in (1.4). 

Let us now consider a general picture of a Lagrangian flow system of equations. Suppose w is 
a contravariant vector quantity attached to a particle following a flow along characteristic paths 
dx/dt = u(x, t) of a velocity field u. Let us consider the abstract Lagrangian flow equation 

Dw . D d _, . „. 

D-t = m + u - v > (L7) 

where the material derivative has its standard definition, and that, in turn, a satisfies the 
Lagrangian equation 

D 2 w Da . , 

-=-5- = = 6(as, t) . 1.8) 

Dt 2 Dt v ' ' K ' 

So far, these are just kinematic rates of change following the characteristics of the velocity 

generating the path x(t) determined from dx/dt = u(x, t). Examples of systems that (1.7) 

might represent are : 

1 . If w represents the vorticity u: = curl u of the incompressible Euler fluid equations then 
a = <jj ■ Vix and divu = 0. With rotation w would be w = u> = Pq 1 ^ + 2fi). 

2. For the barotropic compressible Euler fluid equations (where the pressure p = p(p) is 
density dependent only) then w = u> p = p -1 ^, in which case a = u> p ■ V« and divw = 0. 

3. w could also represent a small vectorial line element S£ that is mixed and stretched by a 
background flow u, in which case a = S£ ■ Vu. For example, following Moffatt's analogy 
with between the magnetic field B in ideal incompressible MHD and vorticity [30] , if w is 
chosen such that w = B, then a = B ■ Vit with div£? = 0. In a more generalized form 
it could also represent the Elsasser variables = u ± B, in which case = ■ Vu 
with two material derivatives. 

4. The semi-geostrophic (SG) model used in atmospheric physics can also be cast in the 
form of (1.7) ; for instance one could choose w = x, a = u and b is computed from the 
SG- model through the semi-geostrophic and a-geostrophic contributions [31, 32, 33]. 

5. For a passive tracer particle with velocity w in a fluid transported by a background velocity 
field u, the particle's acceleration would be a (see [34, 16]). 

In cases (1-3) above if w satisfies the standard Eulerian form 

^ = w-Vu, (1.9) 
then to find b it follows from Ertel's Theorem that [35] 

5fe£W-..v(3s), 

which means that the operators D/Dt and w-V commute for any differentiable function n(x, t). 
Choosing /x = it as in [36], and identifying the flow acceleration as Q(x, t) such that Du/Dt = 
Q(x, t), we have 

D 2 w „ fDu\ 

w-V[ — )=wVQ. (1.11) 



Dt 2 \Dt J 
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In each of the cases (1-3) above Q is readily identifiable and thus we have b 

^ = w-VQ=:b(x,t), (1.12) 

thereby completing the quartet of vectors (it, w, a, b). In §2 it will be shown that knowledge 
of the quartet of vectors (u, w, a, b) determines the quaternion-frame, which is a completely 
natural ortho- normal frame for the Lagrangian dynamics. Modulo a rotation around w, the 
quaternion- frame turns out to be the Frenet-frame attached to lines of constant w. Although 
usually credited to Ertel [35], the result in (1.10), which involves the cancellation of nonlinear 
terms of 0(|w||Vm| 2 ), actually goes much further back in the literature than this; see [36, 37, 
38, 39, 40]. While Ertel's Theorem above enables us to find a b as in cases (1-3), b must be 
determined by other means in case (4). 

1.3 Blow-up in the three-dimensional Euler equations 

The general picture of Lagrangian evolution and the associated quaternion frame is given in §2. 
Thereafter this paper will focus on the three-dimensional incompressible Euler equations (1.1) 
(see §3) and the global existence of solutions (see §4). 

Many generations of mathematicians could testify to the deceptive simplicity of the Euler 
equations. The work of the late Victor Yudovich [41], who proved the existence and uniqueness 
of weak solutions of the two-dimensional Euler equations with uq € L°° on unbounded domains, 
will be remembered as a mile-stone in Euler dynamics. In the three-dimensional case, while many 
special solutions are known in terms of simple functions [16, 17, 18], and powerful results have 
been found on weak and distributional solutions (see Majda & Bertozzi [15]) yet the fundamental 
problem of whether solutions exist for arbitrarily long times or become singular in a finite time 
still remains open. In physical terms, singular behaviour could potentially occur if a vortex is 
resolvable only by length scales decreasing to zero in a finite time. While a review of certain 
aspects of the three-dimensional Euler singularity problem will form part of the later sections of 
this review, the regularity problem for the Navier-Stokes equations will not be considered; the 
interested reader should consult [42, 43, 44]. 

In the first demonstrable case of Euler blow-up, Stuart [45, 46, 47] considered solutions of 
the three-dimensional Euler equations that had linear dependence in two variables x and z; 
the resulting differential equations in the remaining independent variables y and t displayed 
finite time singular behaviour. Stuart then showed how the method of characteristics leads 
to the construction of a complete class of singular solutions [45]. This type of singularity has 
infinite energy because the solution is linearly stretched in the both the x and z directions. 
In a similar fashion, Gibbon, Fokas & Doering [48] considered another class of infinite energy 
solutions whose third component of velocity is linear in z so that the velocity field takes the form 
u = {u\(x,y,t), U2(x,y,t), zj(x,y,t)}. These generalize the Burgers' vortex [18] and represent 
tube and ring-like structures depending on the sign of j(x,y,t). Strong numerical evidence of 
singular behaviour on a periodic x — y cross-section found by Ohkitani and Gibbon [49] was 
confirmed by an analytical proof of blow-up by Constantin [50]. Subsequently Gibbon, Moore 
and Stuart [51] found two explicit singular solutions using the methods outlined in [45]. 

The Beale-Kato-Majda (BKM) theorem [52] has been the main cornerstone of the analysis 
of potential finite energy Euler singularities: one version of this theorem is that ||^||oo must 
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satisfy (see §4 for a more precise statement) 

T 

IMloo dr < oo , (1-13) 

for a global solution to exist up to time T. The most important feature of (1-13) is that it is 
single, simple criterion which is easily monitored. Several refinements of the BKM-Theorem exist 
in addition to those by Ponce [53], who replaced ||^||oo by IjS'lloo, and the BMO-version proved 
by Kozono and Taniuchi [54]. In particular, these take account of the direction in which vorticity 
grows. The work of Constantin [55], and Constantin, Fefferman & Majda [56], reviewed in §4.1, 
deserves specific mention. They were the first to make a precise mathematical formulation of 
how the misalignment of vortex lines might lead to, or prevent, a singularity. §4.2 is devoted 
to a review of the work of Deng, Hou h Yu [57, 58] who have established different criteria on 
vortex lines. In §4.3, quaternions are considered as an alternative way of looking at the direction 
of vorticity [59], which provides us with a different direction of vorticity theorem based on the 
Hessian matrix of the pressure (1.4). Further discussion and references are left to §4. 

1.4 Definition &z properties of quaternions 

In terms of any scalar p and any 3-vector q, the quaternion q = [p, q] is defined as (Gothic fonts 
denote quaternions) 

3 

q = \p, q] = pi - ^2 qi<Ti ' ( L14 ) 
i=i 

where {<ti, a 2 , 03} are the three Pauli spin- matrices defined by 

-00' -(-Si)- ^(°-°)' (L15) 

and which obey the relations a^aj = —SijI — tijk&k ■ I is the 2x2 unit matrix. These relations 
then give a non-commutative multiplication rule 

qi © q2 = [P1P2 - <?i • <?2, PiQ2 + P2Qi + Qi x Q2\ ■ (1-16) 

It can easily be demonstrated that quaternions are associative. One of the main properties of 
quaternions not shared by 3-vectors is the fact that they have an inverse; the inverse of q is 
q* = [p, —q] which means that q © q* = [p 2 + q 2 , 0] = (p 2 + q 2 )[l, 0]; of course, [1, 0] really 
denotes a scalar so if p 2 + q 2 = 1, q is a unit quaternion q. 

A quaternion of the type tt) = [0, w] is called a pure quaternion, with the product between 
two of them expressed as 

tDi © tu 2 = [0, wi] © [0, w 2 ] = [—Wi ■ w 2 , wi x w 2 ] . (1-17) 

In fact there is a quaternionic version of the gradient operator V = [0, V] which, when acting 
upon a pure quaternion u = [0, tt], gives 

V ©u = [-divti, curl it] . (1-18) 

If the field tt is divergence-free, as for an incompressible fluid, then 



L 



V © u = [0, u>] 



(1.19) 
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This pure quaternion incorporating the vorticity will be used freely in future sections. 

It has been mentioned already in Section 1.1 that quaternions are used in the aerospace 
and computer animation industries to avoid difficulties with Euler angles. Here the relation is 
briefly sketched between quaternions and one of the many ways that have been used to describe 
rotating bodies in the rich and long-standing literature of classical mechanics - for more see 
[62]. Whittaker [10] shows how quaternions and the Cayley-Klein parameters [11] are intimately 
related and gives explicit formulae relating these parameters to the Euler angles. Let q = [p, q] 
be a unit quaternion with inverse q* = [p, —q] where p 2 +q 2 = 1. For a pure quaternion r = [0, r] 
there exists a transformation from r — > t' = [0, r'\ 

t' = q©t©q*. (1.20) 

This associative product can explicitly be written as 

r ' = q © r © q* = [0, (p 2 - q 2 )r + 2p(q x r) + 2q(r ■ q)] . (1.21) 
Choosing p = ± cos \9 and q = ±n sin where n is the unit normal to r, we find that 

r' = q © r © q* = [0, r cos 9 + (n x r) sin 9} = 0(9, h)r , (1.22) 

Equation (1.22) is the Euler-Rodrigues formula for the rotation 0(9, h) by an angle 9 of the 
vector r about its unit normal n ; 6 and n are called the Euler parameters. With the choice of 
p and q above q is given by 

q = ±[cos±0, hsm\e}. (1.23) 

The elements of the unit quaternion q are the Cayley-Klein parameters which are related to the 
Euler angles and which form a representation of the Lie group SU(2). All terms in the (1.21) 
are quadratic in p and q, and thus possess the well-known ± equivalence which is an expression 
of the fact that SU(2) covers SO(3) twice. 

To investigate the map (1.20) when p is time-dependent, the Euler-Rodrigues formula in 
(1.22) can be written as 

x'(t)=p®x®p* r = p* ®v'(t) ©p. (1.24) 

Thus t' has a time derivative given by 

x'(t) = p © (p* ©r' ©p) ©p* +p © (p* ©r' ©p) ©p* 

= p © p* © r' + r' © p © p* 

= (p©p*)©r' + r'©(p©p*)* 

= (p©p*)©r'-((p©p*)©r')*, (1.25) 

having used the fact on the last line that because t' is a pure quaternion, r'* = — r'. Because 
P = [Pi q] is °f un h length, and thus pp + qq = 0, this means that p©p* is also a pure quaternion 

£©p* = [0, iOo(t)]. (1-26) 

The 3-vector entry in (1.26) defines the angular frequency fl (t) as fi = 2(— pq + qp — q x q) 
thereby giving the well-known formula for the rotation of a rigid body 



r' = fl x r' . 



(1.27) 
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For a Lagrangian particle, the equivalent of Hq is the Darboux vector T> a in Theorem 1 of §2. 
This Theorem is the main result of this paper and is the equivalent of (1.27) for a Lagrangian 
particle undergoing rotation in flight. 

Finally, it can easily be seen that Hamilton's relation in terms of hyper-complex numbers 
i 2 = j 2 = k 2 = ijk = — 1 will generate the rule in (1.16) if q is written as a 4-vector q = 
p + iqi + jq2 + kq%. Sudbery's paper is still the best source for a study of the functional 
properties of quaternions [60]; he discusses how various results familiar for functions over a 
complex field, such as the Cauchy-Riemann equations, Cauchy's Theorem and integral formula, 
together with the Laurent expansion (but not conformal mappings) have their parallels for 
quaternionic functions. More recent work on further analytical properties can be found in [61]. 

2 Lagrangian evolution equations and an ortho-normal frame 

This section sets up the mathematical foundation concerning the association of quaternion 
frames and can be found in the paper by Gibbon and Holm [62]. Let us repeat the Lagrangian 
evolution equations for a vector field w satisfying (1.7) and (1.8) 

Dw , . Da . . . 

— = a(x,t), — = b(x,t). (2.1) 




Figure 1: The dotted line represents the tracer particle (•) path moving from (xi,ti) to (x%,t-2). 
The solid curves represent lines of constant w to which w is a unit tangent vector. The orientation of 
the quaternion- frame (w, x a , w x x a ) is shown at the two space-time points; note that this is not the 
Frenet-frame corresponding to the particle path but to lines of constant w. 

Given the Lagrangian equations in (2.1), define the scalar a a and the 3-vector x a as5 

a a = \tv\~ 1 (w ■ a) , x a = |*u| _1 (io x a) , (2.2) 

Moreover, let and Xb be defined as in (2.2) for a a and x a with a replaced by b. The 3-vector 
a can be decomposed into parts that are parallel and perpendicular to w (and likewise the same 
for b) 

a = a a w + Xa xw = [a a , x a \ © [°> w \ . ( 2 - 3 ) 

and thus the quaternionic product is summoned in a natural manner. By definition, the growth 
rate a a of the magnitude \w\ obeys 

D\w\ , , . t , 

Dt = aM ' (2 ' 4) 

5 The role of null points w = is not yet clear although, as §3 shows, this problem is neatly avoided by the 
Euler fluid equations. It has been discussed at greater length in [62]. 
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while the unit tangent vector w = ww 1 satisfies 

Dw 



Dt 



X„x-u>. (2.5) 



Now identify the quaternions 6 

q a = [ot a , Xa\ > <\b = [oib, Xb\ , (2-6) 

and let to = [0, w] be the pure quaternion satisfying the Lagrangian evolution equation (2.1) with 
q a defined in (2.6). Then the first equation in (2.1) can automatically be re-written equivalently 
in the quaternion form 

Dto 

— = [0, a] = [0, a a w + Xa x w] = q a © to . (2.7) 

Moreover, if a is differentiable in the Lagrangian sense as in (2.1) then it is clear that a similar 
decomposition for b as that for a in (2.3) gives 

D 2 to 

_ r = [0, b] = [0, a b w + Xb x w] = q b ® to . (2.8) 
Using the associativity property, compatibility of (2.8) and (2.7) implies that (\w\ / 0) 



Dt 

which establishes a Riccati relation between q a and q b 

Dq a 



+ q ®q -q6 ®W = 0, (2.9) 



Dt 



+ qa©qa = q b - (2.10) 



This relation is closely allied to the ortho- normal quaternion-frame 7 (w, % a , w x x a ) whose 
equations of motion are given as follows : 

Theorem 1 [62] The ortho-normal quaternion-frame (w, x a i & x Xa) e SO (3) has Lagrangian 
time derivatives expressed as 

Dw 

= T> ab xw, (2.11) 



Dt 

D(W X Xa) 

Dt 

DXa 



T> ab x(wx Xa) , (2-12) 

VabXXa, (2-13) 



Dt 

where the Darboux angular velocity vector T> ab is defined as 

T>ab = Xa + — W , C b = W ■ {Xa X Xb) ■ ( 2 - U ) 

Xa 

Remark: The analogy with the formula for a rigid body is obvious when compared to (1.27). 
but the Darboux angular velocity vector T> ab is itself a function of x j w an d other variables 
and sits in a two-dimensional plane. In turn this is driven by c b = w ■ (x a x x b ) f° r which b 



6 Dropping the a , b labels and normalizing, the Cayley-Klein parameters are q = [a, xK*^ 2 + X 2 ) ^ 2 - 
7 According to Hanson [2] the the quaternion-frame is similar to the Bishop-frame in computer graphics. 
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must be known. Given this it may then possible to numerically solve equations (2.11) - (2.14) 
for the particle paths. 

Proof : To find an expression for the Lagrangian time derivatives of the components of the 
frame (w, x a , w x x a ) requires the derivative of x a - To find this it is necessary to use the fact 
that the 3-vector b can be expressed in this ortho-normal frame as the linear combination 

w~ 1 b = a b w + c b x a + d b (w xx a )- ( 2 - 15 ) 
where c b is defined in (2.14) and d b = — (x a ' Xb)- The 3-vector product Xb = w _1 (w x b) yields 

X b = c b (w x x a ) ~ d b Xa • (2-16) 

To find the Lagrangian time derivative of x a i we use * ne 3-vector part of the equation for the 
quaternion q a = [a a , x a ] m Theorem 1 

^ = -2a aXa + X b , ^ = -2a aXa -d b , (2.17) 

where Xa = \Xa 

|. Using (2.16) and (2.17) there follows 

^ = C bX -\w X Xa) , D{ ™* t =XaW~ CbXa'Xa , (2.18) 

which gives equations (2.11)-(2.14). ■ 

How to find the rate of change of acceleration represented by the 6-vector is an important 
question regarding computing the paths of passive tracer particles where b is not known through 
Ertel's Theorem. The result that follows describes the evolution of q b in terms of three arbitrary 
scalar s. 

Theorem 2 [62] The Lagrangian time derivative of q b can be expressed as 

= qa©qb + Aiq 6 + A 2 q a + A 3 I, (2.19) 
where the \i(x, t) are arbitrary scalars (I = [1, 0]). 

Proof: To establish (2.19), we differentiate the orthogonality relation Xb ' = an d use the 
Lagrangian derivative of w 

Dxb 

-jyj- = Xa x Xb + s o , where s = ^x a + A Xb ■ ( 2 - 2 0) 

so lies in the plane perpendicular to w in which x a an d Xb a ^ so ne an d A* = ^{x, t) and A = A(aj, t) 
are arbitrary scalars. Explicitly differentiating Xb = w~ 1 (w x b) gives 

w^w (x a -b) + s = -a a x b ~ a b x a + w^w (x a ■ b) + uT 1 x , (2.21) 

which can easily be manipulated into 

w x |^ - a b a - a a b| = w s . (2.22) 

This means that 

Db , 

— = a b a + a a b + s x w + ew , (2.23) 
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where e = e(x, t) is a third unknown scalar in addition to \x and A in (2.20). Thus the Lagrangian 
derivative of a b = w~ x (w ■ b) is 

^ = aa b + Xa -Xb + e. (2.24) 

Lagrangian differential relations have now been found for Xb an d but at the price of intro- 
ducing the triplet of unknown coefficients fi, A, and e which are re-defined as 

A = a a + Ai , n = a b + A 2 , e = -2x a ■ Xb + + Aia fe + A 3 . (2.25) 

The new triplet has been subsumed into (2.19). Then (2.20) and (2.24) can be written in the 
quaternion form (2.19). ■ 

3 Quaternions and the incompressible 3D Euler equations 

The results of the previous section on Lagrangian flows are immediately applicable to the in- 
compressible Euler equations, but to present them in this manner is actually to do so in the 
chronologically reverse order in which they were first developed. Looking ahead in this section, 
the variables a and x in (3-4) for the Euler equations, and the two coupled differential equations 
that they satisfy (3.10), were first written down almost ten years ago in [63, 64] without the 
help of quaternions. It was then discovered in [65] that these equations could be combined to 
form a quaternionic Riccati equation. Finally, the more recent paper [59], in combination with 
[62], put all these results in the form expounded in this present paper. Because data for the 
three-dimensional Euler equations gets very rough very quickly it should be understood that all 
our manipulations are formal. 

In §2 it was shown that a knowledge of the quartet of vectors (u, w, a, b) is necessary to be 
able to use the results of Theorem 1. With w = u> and u = curl it the vortex stretching vector 
is a = u ■ Vw. Thus the w- and it-fields are not independent in this case. Within a = uj ■ Vw, 
the dot-product of u> sees only the symmetric part of the velocity gradient matrix Vtt, which 
is the strain matrix Sij = \ (uij + Uj^) defined in (1.6). With a = u ■ Vtt = Su, the triad of 
vectors is 

(it, (/;. a) = (it, oj, Su)) . (3.1) 

To find the b-field, Ertel's Theorem of §1.2 comes to the rescue. The Du/Dt within the right 
hand side of (1.1) (with w = u>) obeys Euler's equation Du/Dt = —Vp, so we have 

where P is the Hessian of the pressure defined in (1.4). The quartet of vectors necessary to 
make Theorem 1 work is now in place 

(u, w, a, b) = (u, cj, Su>, —Pu>) . (3.3) 

The table below discusses three quartets (u, w, a, b) for the Euler fluid equations : 



u 


7(7 


a 


b 


Material Deriv 


Euler 


X 


u 


-Vp 


(1.2) 


Euler 


U 


-Vp 


? 


(1.2) 


Euler 


u> 


Su 


-Pu 


(1.2) 
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Table 1 : The entries above are three of the possibilities for finding a fo-field given the triplet (u, w, a). 
The third line is the result (3.3) while b is unknown for the second line. 

Using the definitions in §2 the scalar a and the 3-vector \ are defined as 

a = Cj ■ Suj , x = cjxSu>, (3.4) 

together with the definitions for a p and x P 

a p = u ■ Pu , X P = wxPw. (3.5) 

a in (3.4) is now identified as the same as that in Constantin [55] who has expressed it as an 
explicit Biot-Savart formula 8 , a = Su> can be decomposed into parts that are parallel and 
perpendicular to u 

Su> = au> lxxu=[a,x]i) [0, u>] . (3.6) 
By definition, the growth rate a of the scalar magnitude |u;| and the unit tangent vector u obey 

D \ U \ i I D<2 > - /o^ 

— - — = a \uj , — — = yxw, (3.7) 
Dt 1 1 Dt y 1 

which show that a drives the growth or collapse of vorticity and x determines the rate of swing 
of u> around Su. Now identify the quaternions 

q = [a, X], <\p = K>, X P ] ■ (3-8) 

The equivalent of the Riccati equation (2.10) is 9 

^ + q©q + q P = 0, (3.9) 
which, when written explicitly in terms of a~x, becomes 

^ + a 2 - X 2 + a p = 0. R£ + 2a X + X P = 0- (3.10) 

In Theorem 1 we need to use b = —Pu to calculate the path of the ortho-normal quaternion- 
frame (o>, d> x x)- Specifically we must solve 

Du 

— = T>x&, (3.11) 
D <* X *> = X>x(ax*), (3.12) 



Dt 



j£ = Vx X , (3.13) 

where the Darboux angular velocity vector T> is defined as 

T> = x+-&, c p = -u>-( X xx P ), (3.14) 

The pressure Hessian contributes to the angular velocity T> through the scalar coefficient c p . 
To compute the fluid particle paths one would need data on the pressure Hessian P as well as 



8 Everywhere in [55, 56, 66, 67] the unit vector of vorticity is designated as £ whereas here we use d>. 
9 In principle (3.9) can be linearized to a zero-eigenvalue Schrodinger equation in quaternion form with q p as 
the potential, although it is not clear how to proceed from that point. 
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the vorticity oj and the strain matrix S. It is here where the fundamental difference between 
the Euler equations and a passive problem is made explicit. For the Euler equations the 6-field 
containing P is not independent of w = u> but is connected subtly and non-locally through the 
elliptic equation for the pressure (1.5) which we repeat here 

-TrP = Tr{S 2 )-\oj 2 . (3.15) 




Figure 2: The equivalent of Figure 1 but for the Euler equations with the dotted line representing an 
Euler fluid particle (•) path moving from (x\,ti) to (a^,^)- The solid curves represent vortex lines to 
which Cj is a unit tangent vector. The orientation of the quaternion-frame (u>, x, d> x x) is shown at the 
two space-time points; note that this is not the Frenet-frame corresponding to the particle path. 



Theorem 2 expresses the evolution of q p 

Dq p 



Dt 



q ©q p + Aiq p - A 2 q - A 3 I, (3.16) 



in terms of the arbitrary scalars \i(x, t). How these can be determined or handled in terms of 
the incompressibility condition is not clear. 



4 The BKM Theorem & the direction of vorticity 

Three-dimensional Euler data becomes very rough very quickly ; thus, understanding how vor- 
ticity grows and in what direction, are fundamental questions that have yet to be definitively 
answered. Clearly the vortex stretching term cj • Vm = Su>, and the alignment of uj with 
the eigenvectors e.{ of S, play a fundamental role in determining whether or not a singularity 
forms in finite time. Major computational studies of this phenomenon can be found in Bra- 
chet et al. [68, 69]; Pumir k Siggia [70]; Kerr [71, 72]; Grauer et al. [73], Boratav & Pelz [74], 
Pelz [75] , and Hou k, Li [76] . Studies of singularities in the complex time domain of the two- 
dimensional Euler equations can be found in Pauls, Matsumoto, Frisch k Bee [77] where an 
extensive literature is cited. 

The BKM-theorem [52] is the key result in studying the growth of Euler vorticity and possible 
singular behaviour. The domain D C M 3 in Theorem 3 is taken to be a three-dimensional periodic 
domain for present purposes, which guarantees local existence in time of classical solutions (Kato 
[78]), although it is applicable on more general domains than this. One version is {H s denotes 
the standard Sobolcv space): 
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Theorem 3 (Beale, Kato and Majda [52]): On the domain D = [0, L\^ er there exists a global 
solution of the Euler equations, u G C([0, oo]; H s ) n C 1 ([0, oo];£P _1 ) for s > 3 if, for every 
T > 

T 

1 1 w 1 1 £,«> (u) dr < oo . (4.1) 

The result can be stated the opposite way which is that no singularity can form at T without 
Jq ||u>||loo(b) At = oo. Theorem 3 has direct computational consequences. In a hypothetical 
computational experiment if one finds vorticity growth ||u>||l°°(d) ~ (T — t)~ 7 for some 7 > 0, 
then the theorem says that 7 must satisfy 7 > 1 for the observed singular behaviour to be real 
and not an artefact of the numerical computations. The reason is that if 7 is found to lie in 
the range < 7 < 1 then ||ci>||lo°(b) blows up whereas its time integral does not, thus violating 
the theorem. Of the many numerical calculations performed on Euler that by Kerr [71, 72], 
using anti-parallel vortex tubes as initial data, was the first to see 7 pass the threshold with 
a critical value of 7 = 1, followed by Grauer et al. [73], Boratav & Pelz [74] and Pelz [75]. 
Recent numerical calculations by Hou & Li [76], however, have contradicted the existence of a 
singularity: see [79] for a response and a discussion of the issues. To fully settle this question will 
require more refined computations in tandem with analysis to understand the role played by the 
direction of vorticity growth. As indicated in §1, the work of Constantin, Fefferman & Majda 
[56] (see also Constantin [55]) was the first to make a precise mathematical formulation of how 
smooth the direction of vortex lines have to be that might lead to, or prevent, a singularity. §4.1 
is devoted to a short review of this work. Further papers by Cordoba & Fefferman [80], Deng, 
Hou & Yu [57, 58] and Chae [66, 67] are variations on this theme. This approach, pioneered in 
[56], lays the mathematical foundation for the next generation of computational experiments, 
after the manner of Kerr [71, 72, 79] and Hou & Li [76], to check whether a singularity develops. 
§4.2 is devoted to a description of the results in the papers by Deng, Hou & Yu [57, 58] who 
have established different criteria on vortex lines. §4.3 is devoted to an alternative direction of 
vorticity theorem proved in [59] based on the quaternion formulation of this paper. 

References and a more global perspective on the Euler equations can be found in the book 
by Majda and Bertozzi [15]. Shnirelman [81] has constructed very weak solutions which have 
some realistic features but whose kinetic energy monotonically decreases in time and which 
are everywhere discontinuous and unbounded and for its dynamics in the more exotic function 
spaces see the papers by Tadmor [82] and Chae [83, 84, 85]. 

4.1 The work of Constantin, Fefferman & Majda 

The obvious question regarding the BKM-criterion is whether the L°°-norm can be weakened 
to L p for 1 < p < 00. This question was addressed by Constantin [55] who placed further 
assumptions on the local nature of the vorticity and velocity fields. Consider the velocity field 

Ui(t) := sup\u(x, t)\, (4.2) 

X 

and the L^-norm of u defined by 

IM|i>c = £- 3 sup / \u(x + y)\d 3 y, (4.3) 

x J\y\<L 



L 
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where L is some outer length scale in the Euler flow which could be taken to be unity. Now 
assume that the unit vector of vorticity is Lipschitz 

\&(x,t)-&(y,t)\<\^£ (4.4) 

for \x — y\ < L and for some length po(t). Then the following result is stated in Constantin [55] 
and re-stated and proved in Constantin, Fefferman & Majda [56]: 

Theorem 4 (Constantin [55]; Constantin, Fefferman & Majda [56]): Assume that the initial 
vorticity cj is smooth and compactly supported and assume that a solution of the Euler equations 
satisfies 

<*>(• , s)|| Moc ( — ] ds <oo, / — — ds <oo. (4.5) 

\Po(s)J Jo Po(s) 

Then 

sup rr-7 .>.„ < oo . (4.6) 

0<t<T ||^(- , t)\\l,loc 

Clearly, if U\ = ||«||oo < 00 ana ll^(' > *)||i, ioc < oo on [0, T] and po is bounded away from zero 
then the BKM theorem says that no singularities can arise. The Lipschitz condition (4.4) can 
be re-expressed to account for anti-parallel vortex tubes [55]. 

Constantin, Fefferman & Majda [56] then considered in more detail how to define the concept 
of "smoothly directed" for trajectories. Consider the three-dimensional Euler equations with 
smooth localized initial data; assume the solution is smooth on < t < T. The velocity field 
defines particle trajectories X(xq, t) that satisfy 

DX 

— = u(X,t) X(x ,0)=x . (4.7) 

The image Wt of a set Wo is given by Wt = X(t, Wo). Then the set Wq is said to be smoothly 
directed if there exists a length p > and a ball < r < \p such that the following 3 conditions 
are satisfied: 

1. For every x <G W * where W * = {x G W * ; |cj (^o)| 7^ 0} and all t <G [0, T), the function 
(*>(■ , t) has a Lipschitz extension to the ball of radius 4p centred at X(xq, t) and 



M=lim sup f ||Vu>(-, t)\\ 2 Loo(Bl )dt 



< oo . (4.8) 



This assumption ensures the direction of vorticity is well-behaved in the neighbourhood 
of a set of trajectories. 

2. The condition 

sup \u)(x, t)\<m sup |u>(x, t)\ (4.9) 

B» r {Wt) B r (Wt) 

holds for all t € [0, T) with m = const > 0. This simply means that this chosen neigh- 
bourhood captures large & growing vorticity but not so much that it overlaps with another 
similar region. 

3. The velocity field in the ball of radius 4/? satisfies 



sup \u(x, t)\ < U(t) := sup \u(x, t)\ < oo , (4.10) 

B ir (W t ) x 



for all t G [0, T). 
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Theorem 5 ( Constantin, Fefferman & Majda [56]) Assume that Wo is smoothly directed as in 
(i)-(iii) above. Then there exists a time r > and a constant V such that 

sup \u>(x,t)\<T sup \u>(x, to)\ (4-11) 

B r (Wt) B p (W t ) 
holds for any < to < T and < t — to < r. 

Condition (ii) may have implications for how the natural length p scales with time as the flow 
develops [72] but more work needs to be done to understand its implications. Cordoba & 
Fefferman [80] have weakened condition (iii) in the case of vortex tubes to 

f U(s)ds= [ \\u(-, s) lUds < oo. (4.12) 
Jo Jo 

4.2 The work of Deng, Hou & Yu 

Deng, Hou & Yu [57] have re-worked probably the most important of the "smoothly directed 
criteria", namely (4.8), from local control over J* ||Vd>(- , t)\\ 2 Loo dt in < t < T to a condition 
on the arc-length s between two points s\ and S2- The first of their two results is : 

Theorem 6 (Deng Hou & Yu [57]): Let x(t) be a family of points such that \us(x(t), t)\ > 
£l(t) = ||u;||oo- Assume that for all t £ [0, T] there is another point y(t) on the same vortex 
line as x(t) such that the unit vector of vorticity u>(x, t) along the line between x{t) and y{t) is 
well-defined. If we further assume that 

divu(s, t)ds < C(T) (4.13) 
\u (x (t), t)\dt < oo , (4.14) 



Sl 



together with 

f-T 



J 

Jo 



10 

then there will be no blow-up up to time T. Moreover, 



e~ C < P^4 <e c . (4.15) 
" Wv{t), t)\~ 

Inequality (4.13) is based on the simple fact that 

d\ijj\ 

= diva; = IcjI diva) + a> • Vlwl = \u> \ diva; -| — (4-16) 

ds 

where u ■ V = ^ is the arc-length derivative. 

The second and more important of the results of Deng, Hou & Yu [58] is based on considering 

a family of vortex line segments Lt along which the maximum vorticity is comparable to the 
maximum vorticity f2(t). Denote by L(t) the arc length of Lt, h the unit normal vector, and k 
the curvature of the vortex line. Furthermore, they define 

Ua(t)= max \(u-Cj)(x,t)-{u-u){y,t)\ , (4.17) 

x,y€L t 

U n (t) = max \u ■ h\ , (4-18) 

L t 

and 

M(t) = max (J|V • u>\\ Loo(Lt) , ||k|| L oc (L() ) . (4.19) 
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Theorem 7 (Deng Hon & Yu [57]) : Let A,B € (0, 1) with B < 1 - A, and C be a positive 
constant. If 

1. uut) + u n (t)<(T-tr A , 

2. M(t)L(t) < C , 

3. L(t)>(T-t) B , 

then there will be no blow-up up to time T . 

In a further related paper Deng, Hou & Yu [58] have changed the inequality A + B < 1 to 
equality A + B = 1 subject to a further weak condition. They also derived some improved 
geometric scaling conditions which can be applied to the scenario when the velocity blows up 
at the same time as vorticity and the rate of blow-up of velocity is proportional to the square 
root of vorticity. This is the worst possible blow-up scenario for velocity field due to Kelvin's 
circulation theorem. 



4.3 The non-constancy of a p & x P - quaternions & the direction of vorticity 

The key relation in the quaternionic formulation of the Euler equations is the Riccati equation 
(3.9) for q = [a(x,t), x(x, t)]. In terms of a and X this gives four equations 

^ = X 2 ~o?-a p , ^f = -2a X - Xp . (4.20) 

Although apparently a simple set of differential equations driven by q p = [a p , x P ]> it is clear 
that q p is not independent of the solution because of the pressure constraint —TrP = it^fcUfc^. 
In consequence it is tempting to think of q p as behaving in a constant fashion. This may be 
true for large regions of an Euler flow but it is certainly not true in the most intense vortical 
regions where vortex lines have their greatest curvature; in these regions the signs of a p and of 
the components of x P may change dramatically [64]. It is because of these potentially violent 
changes that q p could be considered as a candidate for a further conditional direction of vorticity 
theorem along the lines of those in §4.1 and §4.2. Other work where constraints on P appear is 
the paper by Chae [67]. 

The work in [56, 57, 58] shows that Vw needs to be controlled in some fashion in local areas 
where vortex lines have high curvature. In terms of the number of derivatives the Hessian P is 
on the same level and it is in terms of P and the variables associated with it (a p and Xp) that we 
look for control of Euler solutions. From their definitions, it is easily shown that a 2 + x 2 = 1 5"'^' 1 2 
and thus on vortex lines, a = a(X(t, xq), t), (4.20) becomes 

^\Su\ 2 = -a\Su\ 2 + aa p + x-X P (4.21) 

Thus on integration 

\Su>(X(t), t)\ 2 = -2 f T e fo dt'- ft a(-,t')dt> (aa p + x-X P ) W, r) dr . (4.22) 

There are now two alternatives. The first is to make one application of a Cauchy-Schwarz 
inequality and use the fact that a p + x p = \P&\ 2 

\Su>{X{t, x ), t)\ < 2 f T efo «(•-*') dt '~ So a(-,f)cif T )| dr . (4.23) 

Jo 
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This is similar to Chae's result (his Theorem 5.1 in [67]) which is based on control of the time 
integral of \\Su> ■ Pu\\oo, which is derivable from (3.2). 

The second raises an interesting case respecting the direction of vorticity using x P an d 
can be viewed as an alternative way of looking at the direction of vorticity after [56, 57, 58]. 
X p = & x Pu contains u not a; & is thus concerned with the direction of u> rather than its 
magnitude. Firstly we use the fact that |cj| cannot blow-up for a < because D\uj\/Dt = a\u>\ ; 
thus our concern is with a > 0. In the case when the angle between u and Pu> is not zero 

\Su(X(t, xq), t)\ < 2 / \ Xp (; r)\dr. (4.24) 
Jo 

If the right hand side is bounded then Euler cannot blow up, excepting the possibility that 
\PCj\ blows up simultaneously as the angle between a> and Ptb approaches zero while keeping 
X p finite; under these circumstances J * IxJdr < oo, whereas Jg |a p |dr — > oo and thus blow-up 
is still theoretically possible in that case. The result does not imply that blow-up occurs when 
collinearity does; it simply implies that under condition (4.24) it is the only situation when it 
can happen. Ohkitani [36] and Ohkitani and Kishiba [40] have noted the collinearity mentioned 
above; they observed in Euler computations that at maximum points of enstrophy, u> tends to 
align with the eigenvector corresponding to the most negative eigenvalue of P. Expressed over 
the whole periodic volume we have : 



Theorem 8 (Gibbon, Holm, Kerr & Roulstone [59]) : On the domain B = [0, L] z per there exists 
a global solution of the Euler equations, u G C([0, oo];H s ) n C 1 ([0, oo]; i? s_1 ) for s>3if for 
every T > 

[ ||Xp(-,r)|| 00 dT<oo. (4.25) 
J o 

excepting the case where u becomes collinear with an eigenvector of P at T. 

5 A final example: the equations of incompressible ideal MHD 

The Lagrangian formulation of §2 can be applied to many situations, such as the stretching of 
fluid line-elements, incompressible motion of Euler fluids and ideal MHD (Majda & Bertozzi 
[15]). We choose ideal MHD in Elsasser variable form as a final example; another approach to 
this can be found in [86] . The equations for the fluid and the magnetic field B are 

^ = B-VB-Vp, (5.1) 

DB , s 

— = B-Vu, (5.2) 

together with divit = and div£? = 0. The pressure p in (5.1) is p = pf + \B 2 where pj is the 
fluid pressure. Elsasser variables are defined by the combination [30] 

v ± = u±B. (5.3) 

The existence of two velocities v ± means that there are two material derivatives 



23/01/07 



2nd version for Russian Math Surveys 



19 



In terms of these, (5.1) and (5.2) can be rewritten as 

D ± v* 



Dt 



= -Vp, 



with the magnetic field B satisfying (divv^ = 0) 

D ± B 



Dt 



B ■ Vv ± . 



(5.5) 



(5.6) 



Thus we have a pair of triads (v ± , B, a^) with a ± = B ■ Vv ± , based on Moffatt's identification 
of the .B-field as the important stretching element [30]. From [65, 59] we also have 



D ± gT 
Dt 



PB 



(5.7) 



where b^ 1 = —PB. With two quartets (t» ± , B, ,b), the results of §2 follow, with two La- 
grangian derivatives and two Riccati equations 



(5.8) 



In consequence, MHD-quaternion-frame dynamics needs to be interpreted in terms of two sets of 
ortho-normal frames (^B, % ± , B x acted on by their opposite Lagrangian time derivatives. 



D^B 
Dt 



£>t x B 



Dt A ' 

where the pair of Elsasser Darboux vectors T> T are defined as 



6 Conclusion 



X 



x- 



4 = B ■ [x ± x ( XpB + « ± X T )] • 



(5.9) 
(5.10) 
(5.11) 

(5.12) 



The well-established use of quaternions by the aero/astro-nautics and computer animation 
communities in the spirit intended by Hamilton gives us confidence that they are applicable 
to the 'flight' of Lagrangian particles in both passive tracer particle flows and, in particular, 
three-dimensional Euler flows. An equivalent formulation for the compressible Euler equations 
([46, 47]) may give a clue to the nature of the incompressible limit [87]. The case of the barotropic 
compressible Euler equations and other examples are given in the summary below in Table 2 : 



System 


u 


w 


a 


b 


Material Dcriv 


incompressible Euler 


u 


X 


u 


-Vp 


D/Dt 


incompressible Euler 


u 


u 


Su 


-Pu 


D/Dt 


barotropic Euler 


u 


u/p 


(w/p) • Vit 


-(u j /p)d j (pd i p) 


D/Dt 


MHD 


v ± 


B 


B ■ VvT 


-PB 


D ± /Dt 


Mixing 


u 


S£ 


Si ■ Vu 


-Pd£ 


D/Dt 
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Table 2 : The entries display various examples of the use of Ertel's Theorem in closing the quartet of 
vectors (u, w, a, b). For ideal MHD, D ± /Dt is defined in (5.4). 

Whenever quaternions appear in a natural manner, it usually a signal that the system has 
inherent geometric properties. For the Euler equations, it is significant that this entails the 
growth rate a and swing rate x of the vorticity vector, the latter being very sensitive to the 
direction of vorticity with respect to eigenvectors of S. To elaborate further, consider a Burgers' 
vortex which represents a vortex tube [18]. An eigenvector of S lies in the direction of the 
tube-axis parallel to u in which case x = & x = 0. However, if a tube comes into close 
proximity with another then they will bend and maybe tangle. As soon as the tube-curvature 
becomes non-zero along a certain line-length then along that length. Likewise this will 

also be true for vortex sheets that bend or roll-up when in close proximity to another sheet. 
The 3- vector x 1S therefore sensitively and locally dependent on the vortical topology. In fact at 
each point its evolution is most elegantly expressed through its associated quaternion q, which 
must satisfy (see (3.9)) 

^ + q©q + q P = 0. (6.1) 

To fully appreciate the power of the method the pressure field must necessarily appear explicitly 
in the form of its Hessian through q p although this runs counter to conventional practice in fluid 
dynamics where it is usually removed using Leray's projector. The pressure Hessian appears in 
the material derivative of the vortex stretching term, through the use of Ertel's Theorem, as 
the price to be paid for cancelling nonlinearity 0(\u\ |Vti| 2 ). In fact, the effect of the pressure 
Hessian on the vorticity stretching term is subtle and non-local. Therefore, while it is tempting 
to discount the pressure because it disappears overtly in the equation for the vorticity, covertly 
it may arguably be one of the most important terms in inviscid fluid dynamics. 

There are, of course, stationary solutions of (6.1) one of which is x = Xp = with a = ao and 
The Burgers' vortex is a solution of this type: see [64, 65]. Having laid much stress 
in §4.3 on the non-constancy of a p and x P m intense, potentially singular regions, nevertheless 
let us to try to determine the simplest generic behaviour of a and x from (4.20) when a p and x P 
are constant; for example, a near-Burgers' vortex. To do this let us consider the four equations 
which come out of (6.1), as in (4.20), and think of them as ordinary differential equations on 
particle paths X(t,Xo) 

a = X 2 ~ o> 2 - a p , x = -2ax - C p . (6.2) 

In regions of the a — x phase plane where a p = const, C p = x " X P = const there are 2 critical 
points 

(a, X ) = (±ao, xo) 2a 2 = a p + [a 2 p + C^ 2 (6.3) 

Thus there are two fixed points; one with ao > (stretching), which is a stable spiral, and one 
with ao < (compression); both have a small and equal value of xo- The point with ao < is 
an unstable spiral while ao > is stable. Perhaps it is a surprise that it is the stretching case 
that is the attracting point although it should also be noted that these equations without the 
Hessian terms have arisen in Navier-Stokes turbulence modelling [88]. 

Finally, the existence of the relation (6.1), and its more general Lagrangian equivalent (3.9), 
is the key step in proving Theorem 1, from which the frame dynamics is derived. Moreover, for 
the three-dimensional Euler equations, (6.1) is also the key step in the proof of Theorem 8. 
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